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Abstract 



We develop a new method in lattice QCD to calculate the form factor 
J^B^D*(1) at zero recoil. This is the main theoretical ingredient needed to de- 
termine \Vcb\ from the exclusive decay B D*W. We introduce three ratios, 
in which most of statistical and systematic error cancels, making a precise 
calculation possible. We fit the heavy-quark mass dependence directly, and 
extract the I/tuq and three of the four I/tuq corrections in the heavy-quark 
expansion. In this paper we show how the method works in the quenched 
approximation, obtaining ^^^^.(l) = O.giSt":!]?^ ± o.Oie^nj^+O-OOO+o.o^^ 

where the uncertainties come, respectively, from statistics and fitting, match- 
ing lattice gauge theory to QCD, lattice spacing dependence, light quark mass 
effects, and the quenched approximation. We also discuss how to reduce these 
uncertainties and, thus, to obtain a model-independent determination of \ Vcb\- 

PACS numbers: 12.38. Gc, 12.15.Hh, 13.20.He 
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I. INTRODUCTION 



In flavor physics the Cabbibo-Kobayashi-Maskawa (CKM) matrix element Vcb plays an 
important role. Much of the phenomenology of CP violation centers around the unitarity 
triangle, and a precise value of \Vcb\ is needed to locate the triangle's apex in the complex 
plane. As a fundamental parameter of the Standard Model, Vet sometimes appears in unex- 
pected places. For example, the Standard Model prediction of the K^-K^ mixing parameter 
ex is very sensitive to \Vch\ 

The determination of is made through inclusive and exclusive semileptonic B decays, 
but at present both methods are limited by theoretical uncertainties. The inclusive method 
requires a reliable calculation of the total semileptonic decay rate of the B meson, which 
can be done using the heavy quark expansion Ultimately this method is limited by 

the breakdown of local quark-hadron duality, which is difficult to estimate. The exclusive 
method, on the other hand, requires a theoretical calculation of the form factor Tb^d* of 
B D*lu decay. In this paper we take a step towards reducing the uncertainty in the 
exclusive method, by devising a precise method to compute the form factor at zero recoil in 
lattice QCD. 

The differential rate for the semileptonic decay B —>■ D*lui is given by 

^ = ^ml,{mB-mo*f^Ji;^^g{w)\VM^B^D*{w)\\ (1.1) 
aw Ati'^ 

where w = v' ■ v is the velocity transfer from the initial state (with velocity v) to the final 
state (with velocity v'). The velocity transfer is related to the momentum q transferred 
to the leptons by = m\ — 2wmBm£)* + m|,,, and it lies in the range 1 < w < {m% + 
mjj*) / 2-171 Bfn D* ■ The function 

w + 1 / 8w{w - l)mBmD* \ 

G{w) = —— \5w + 1 + — 1.2 

12 V [mB-mD*y J 

has a kinematic origin, with ^(1) = 1. Thus, given the form factor J^b^d*{w)., 
the measured decay rate to determine \Vcb\- 



One makes use of the zero-recoil point w = 1, even though the phase-space factor a/w^ — 1 
suppresses the event rate, because then theoretical uncertainties are under better control. 
For w > 1, J^B-^D*{w) is a linear combination of several form factors of S ^ D* transitions 
mediated by the vector and axial vector currents. At zero recoil, however, 

J-B^Z)*(1) =/iA,(l), (1.3) 
where /i^i is a form factor of the axial vector current namely. 



{D*{v)\A^\B{v)) = i^2mB 2mD' e'^/iAi(l). (1.4) 

More importantly, heavy-quark symmetry plays an essential role in constraining /i^^(l), 
leading to the simple heavy quark expansion 



/lAi(l) = r]A 



1 _ + 2£ 



i2mcY 2mc2mb {2mbY 



:i.5) 
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including all terms of order l/rriQ. In Eq. ( p..5|) , r]A is a short- distance radiative correction, 
which is known at the two-loop level , and the is are long-distance matrix elements of the 
heavy-quark effective theory (HQET).Q Heavy-quark symmetry normalizes the leading term 
inside the bracket to unity |§| and, moreover, forbids terms of order l/mg The l/rriQ 
corrections are formally small — (A/2mc)^ ~ 4% — but one would like to reach better precision 
on \Vcb\, so these terms cannot be neglected. 

There have been mainly two different methods used to estimate the l/rriQ terms in 
Eq. (|1.5|) , but neither has achieved a model independent calculation. One involves using 
a quark model |^,T^] to estimate the is. The other employs the zero- recoil sum rule [IT . 
Although based on a rigorous upper bound [1^, to make a prediction of this 



approach requires an assumption on the effects of higher excited states in the sum rule. 
Thus — ^just as with quark models — it is difficult to estimate, let alone reduce, the uncertainty 
associated with the estimate. 

In this paper we take a step towards reducing the theoretical uncertainty by using lattice 
QCD to calculate hA^{l) = J-'b^d*{^)- Lattice QCD is, in principle, model independent, 
although here we work in the quenched approximation. The quenched approximation is not 
less rigorous than the methods used in Refs. [|10|,|lT|. From our point of view, however, the 



main advantage of the quenched approximation is that it allows us to learn how to control 
and estimate all other lattice uncertainties. With a proven technique, it is conceptually 
straightforward, if computationally demanding, to carry out a calculation in full QCD. 

Until now three obstacles prevented even quenched lattice calculations of /iAi(l) fo the 
needed precision. First, a direct Monte Carlo calculation of the matrix element in Eq. ( |1.4|) 
suffers from a statistical error that is too large to be interesting. Second, the normaliza- 
tion of the lattice axial vector current was uncertain, being limited by a poorly converging 
perturbation series. Finally, early works ||13[ used ad hoc methods for heavy quarks on the 
lattice, which entailed a poorly controlled extrapolation in the heavy quark mass. We have 
devised methods to circumvent all three obstacles. The first two are handled with certain 
double ratios of correlation functions, in which the bulk of statistical and systematic un- 
certainties cancel |14]. The third obstacle — the problem of heavy-quark lattice artifacts — is 
overcome by using a systematic method for treating heavy quarks on the lattice, based on 
Wilson fermions |15|. This obstacle could also be overcome using lattice NRQCD [T^, as in 



the work of Hein et at . 

In our work on the form factor in the decay B 

role was played by the double ratio of matrix elements 



7^^ 



{D\c-f%\B){B\b-f^c\D) 
{D\c-i^c\D){B\h%\B) 



|/^+(l)P 



Dlv at zero recoil, a central 



:i.6) 



where 



{D{v)\V^'\B{v)) =iV2mB 2mD v^'h+il] 



[1.7) 



^In the HQET literature, the is are often called "hadronic parameters" , because they are viewed 
as incalculable. In a QCD context, however, the are not free parameters, but calculable matrix 
elements. 
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In Ref. [|T4[ we studied the heavy-quark mass dependence of using a fit to obtain the 

l/rriQ and I/ttlq corrections. In this work we employ this double ratio and two similar ones. 
The first additional double ratio is 



7^l 



{D*\c-f%\B*){B*\Irf^c\D*) 
{D*\c'j^c\D*){B*\b-f%\B*) 



Ml] 



(1.8) 



where the pseudoscalar mesons B and D, and their form factor /i+(l), are replaced with the 
vector mesons B* and D*, and their form factor hi{l): 



{D*{v)\V''\B*{v)) = W2mB* 2mD* e' ■ ev^hi{l] 



The second additional double ratio is 



7^. 



{D*\cj,j,b\B){B*\h,j5c\D) 
{D*\c^,^,c\D){B*\hjl5b\B) 



(1) 



l,D~ 



:i.9) 



;i.io) 



leading to a 
the double ratios overcome 



where the axial vector current mediates pseudoscalar-to-vector transitions, 

double ratio of the form factor h^i- As stressed in Ref. |^ ^ 

two of the obstacles in the lattice calculation, because numerator and denominator are so 
similar. Statistical fiuctuations in the numerator and denominator are very highly correlated 
and largely cancel in the ratio. Also, most of the normalization uncertainty in the lattice 
currents cancels, leaving only a residual normalization factor that can be computed reliably 
in perturbation theory |l^. Indeed, all uncertainties scale as TZ — 1, rather than as TZ. 

Note that the double ratio TZa^ does not yield the desired form factor h^'^^* , but instead 
the combination Ka^ , which is itself a double ratio of form factors. One can, however, extract 
hAi{l) from the three double ratios TZi, and TZa^, at least to the order in the heavy- 
quark expansion given in Eq. ( p.. 51 ). This possibility follows from the heavy quark expansions 
for h+{l) and h^{l) §§, 



^+(1) = Vv 
hi{l) = r]v 



1 - 



■V 



— -—)' 



:i.ii) 
:i.i2) 



IS 



and comparing to Eq. ( [1.51 ). In /i+(l) and hi{l) the absence of terms of order l/mq 
easily understood, because charge conservation requires /i+(l) = hi{l) = 1 when rric = rrih, 
and because the matrix elements defining them are symmetric under the interchange rric 
rrih. Similarly, the heavy-quark expansion of the form factor ratio /i^j(l), obtained from 
TZai, is 



/iAi(l) = fjA 



1 - 



'A 



1 



2mr 



2mhJ 



(1.13) 



which follows immediately from Eq. ( |1.5| ), defining fj\ = riAcbriAb<: / rjA'^'^fjAbb . Hence, by varying 
the heavy quark masses in the lattice calculation of the double ratios 7Z+, TZi, and TZa^, one 
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can extract £p, and respectively. Then, /iAi(l) = ^b^d*{X) can be reconstituted 
through Eq. ([L5|) . 

A key to this method is that heavy-quark symmetry requires the quantities ^p and to 
appear in Eq. (|1.5|) , as well as in Eqs. ( p..llj ) and ( |1.12|) 0,^. A simple argument explains 



why. For each form factor there are three possible terms at order l/mQ — ^/^l-, ^/^ti and 
l/iricmb — and each multiplies an HQET matrix element. For and hi{l) the particular 

form of the expansions is restricted by the b ^ c interchange symmetry, so only one HQET 
matrix element can appear in each case: ip for and iy for hi{l). Interchange symmetry 

does not apply to the B — > D* transition, however, so three HQET matrix elements are 
needed in the expansion of /^^^^(l), Eq. (|1.5| ). Two of them, however, coincide with ip 
and £v- If one flips the spin of the charmed quark in the B ^ D transition in Eq. (|1 . 7| ) , one 
obtains the B ^ D* transition in Eq. ( |1.4| ), and in the limit of infinite charmed quark mass 
the matrix elements are identical, by heavy-quark spin symmetry. Consequently, the 
term in Eq. ( |1.5| ) must be the same as that in Eq. ( |1 . 1 1| ) , namely £p/(2mf,)^. The same logic 
applied to the b quark's spin, starting from the B* D* transition in Eq. ( |1.9|) , implies 
that the term in Eqs. ( |1.5| ) and ( |1.12| ) must be the same, namely iv/{'2mc^'^ 



At order I/itlq there are, in general, four terms for each form factor. In Sec. we show 
how the same kind of reasoning can be used to extract three of the four terms from the 
l/rriQ behavior of the three double ratios. Including these corrections not only reduces the 
systematic error of the heavy quark expansion, but also reduces our statistical error, because 
fitted values for the quadratic and cubic terms are correlated. 

In the remainder of this paper we describe the details of our lattice calculation of 
JFb^£)*(1) = hAi{l), as sketched above. Discretization effects are studied by repeating 
the analysis at three different lattice spacings. The dependence on the light quark mass 
is expected to be small, which we are able to verify. After a thorough investigation of 
systematic uncertainties, we obtain 



T /I N _ n Qi Q+U.U24 _|_ r, (-+u.uu3+u.uuu+u.uub (^^A\ 

J-B^D*[i^) - U.yid_o.oi7 ± U.Uib_o.oi4-0.016-0.014 {^■^'^) 

where the uncertainties come, respectively, from statistics and fitting, matching lattice gauge 
theory and HQET to QCD, lattice spacing dependence, light quark mass effects, and the 
quenched approximation. A preliminary report of this calculation based on our coarsest 
lattice appeared in Ref. H, reporting J^b^d*{1) = 0.935 ± 0.022t[!:!]i. The change comes 
mostly from the results on two finer lattices, partly from some secondary changes in the 
analysis, and partly from the inclusion of some contributions of order I/itlq. Clearly, these 
central values are indistinguishable within the error bars. 

The paper is organized as follows. In Sec. || we discuss how to combine heavy-quark 
theory and lattice gauge theory to calculate the needed matrix elements; in particular, we 
review how we are able to extract the 1 / tuq corrections |2^ . Section |I| is fairly general and 
much of it also applies to lattice NRQCD. Specific details of our numerical work are given 
in Sec. |ITT|, including input parameters and the basic outputs. The "Fermilab" method for 
heavy quarks requires matching the short- distance behavior of lattice gauge theory to 



QCD, which is discussed in Sec. 0. Section |V| shows a key feature of our analysis, namely 
the direct fitting of the heavy-quark mass dependence to obtain the power corrections in 
Eq. (|1.5[). A detailed discussion of the systematic uncertainties is in Sec. |VI|. Our result. 
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Eq. (|1.14| ), is compared to other methods in Sec. [VI ]| . Section |VIII| contains some concluding 
remarks. 



II. CONTINUUM AND LATTICE MATRIX ELEMENTS 

In this section we discuss how to obtain continuum-QCD, heavy-quark observables from 
lattice gauge theory. Discretization effects of the heavy quarks are a special concern, so 
they are discussed in detail in this section. For the light spectator quark we use well-known 
methods, and we provide details in Sec. p| . 

Discretization effects of the heavy quarks can be controlled by matching the lattice theory 



to HQET [Q. This is possible whether one discretizes the NRQCD effective Lagrangian |]I6 
or one employs the non-relativistic interpretation of Wilson fermions In either case, 

on-shell lattice matrix elements can be described by a version of (continuum) HQET, with 
effective Lagrangian (in the rest frame) 

Ckqet = mih^hy ^ \ \ , (2.1) 

21712 ZniQ 

where hy is the heavy-quark field of HQET, and B is the chromomagnetic field. The 
"masses" mi, m2, and are short-distance coefficients; they depend on the bare couplings 
of the lattice action, including the gauge coupling. Matrix elements are completely indepen- 
dent of mi ^y], so the important coefficients are m2 and ms- The lattice NRQCD action 



has bare parameters that correspond directly to m2 and mg. With Wilson fermions one 
must use the Sheikholeslami-Wohlert (SW) action and adjust mo and csw to tune m2 



and mg. In practice, we tune m2 non-perturbatively, using the heavy-light and quarkonium 
spectra, and mg with the estimate of tadpole- improved, tree-level perturbation theory ||22|| . 
There are also terms of order l/mg in the effective Lagrangian £hqet, but they do not 
influence the double ratios, as discussed further below. 

In this paper we use lattice currents that are constructed as in Ref. ]T5|]. (An analogous 
set of currents can be constructed for lattice NRQCD |^.) We distinguish the lattice 
currents V'^ and from their continuum counterparts V and A'^. We define 



V^' = JZyccZybt^cil^^b (2.2) 



= JZvccZy,,^,iYl5^b (2.3) 



where the rotated field 



= [1 + ad,j ■ £)lat] i^g, (2.4) 

and tpq is the lattice quark field (g = c, b) in the SW action. Here Z^iat is the sym- 
metric, nearest-neighbor, covariant difference operator. In Eqs. ( ^.2[ ) and ( |2.3| ) the factors 
Zyii, q = c,b, normalize the flavor-conserving vector currents. Because for massive quarks 
only Zy can be computed non-perturbatively, we choose to put Zy into the definition of 
the axial current A^. In the work reported in this paper, we do not need to compute the 
factor ^/ZYccZYbb, because it cancels in the double ratios. 



6 



Matching the current to HQET requires further short- distance coefficients: 



T/M^r'-S,,-, B'^%p^t^^K B^fi^fli^p^ 

V -C^,t;cA +■■■, (2.5) 

where the symbol = imphes equahty of matrix elements, and hy and are HQET fields for 
the bottom and charmed quarks. At the tree level the short-distance coefficients C|^*, C]^|, 
and -Bjfj all equal one. The free parameter di in Eq. ( p. 41) can be adjusted to tune l/mag 
to l/mg. In the present calculations, we adjust di with the estimate of tadpole- improved. 



tree-level perturbation theory, as explained in Ref. [|T5[. Further dimension-four operators, 
whose coefficients vanish at the tree level, are omitted from the right-hand sides of Eqs. ( p.5|) 
and (|2]|); they are listed in Ref. [0. 

The description in Eqs. (|2.5| ) and (|2.6|) is in complete analogy with that for the continuum 
currents, namely, 

V -Cv,v cyby — — + ■ ■ ■ , (2.7) 
^ = ^^.^^^^^^^'^ + + ■ ■ • • (2-8) 



The radiative corrections to the short- distance coefficients in Eqs. ( |2.5| ) and ( |2.6| ) differ from 
those in Eqs. (|2.7| ) and (|2.8| ), because the lattice modifies the physics at short distances. On 
the other hand, the HQET operators are the same throughout. 

There are also terms of order l/mq in the effective currents on the right-hand sides of 
Eqs. ( p.5D -( P^ , although for brevity they are not written out. The most important operator 
in each case is 



yl^ ^ CvP±v^P^K ^2 9) 

A, c^P^nli.P^K_ (2.10) 



As the notation suggests, both these currents are correctly normalized at the tree level when 
di is adjusted so that mag = mq, as above. In addition to these l/nicmb currents, there are 
currents of order and Although the latter contribute to the individual matrix 

elements {D^*^\J'^\B^*'>), their contributions drop out of the double ratios. 

In the foregoing discussion, most corrections of order 1 / mg have been handled only in a 
cursory way. Since we aim for the ^/ttlq corrections to the double ratios we must, however, 
discuss how these contributions are incorporated, when the lattice action and currents are 
constructed and normalized along the lines given above. The HQET description of matrix 
elements reveals several sources of such contributions 

1. double insertions of the l/mg terms in the effective Lagrangian Chqet] 
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2. single insertions of the l/mg terms in the effective Lagrangian into matrix elements 
of the l/rriQ terms in the effective HQET currents; 

3. single insertions of genuine l/mg terms in the effective Lagrangian; 

4. matrix elements of genuine l/mg terms in the effective HQET currents. 

The first set of contributions is correctly normalized at the same level of accuracy as the 
l/mg terms of the action. The second set makes no contribution to zero recoil matrix 
elements whatsoever [Q. The third set also makes no contribution at zero recoil, because 
the leading terms in Eqs. ( p^.5D and ( p^.6D are Noether currents of the heavy-quark symmetries 
and, as in the proof of Luke's theorem, first corrections to Noether currents vanish [p5| , p0| 



One is left with the last set, which does contribute to the matrix elements defining the 
form factors. The HQET matrix elements of all dimension-five currents can be reduced to 
Ai and A2, which appear in the heavy-quark expansion of the mass @]. In the double ratios, 
however, the following cancellation (schematically) takes place [pO|j : 

[1 - XiXb/ml - + Xjml)]^ -i_x(J--l-)\ (2.11) 



[l-A(2X,-l)/m2][l-A(2Xb-l)/m2] Vm, 

where A is proportional to Ai or A2, and Xq/ttlq indicates incorrect normalization, while 
l/mgrrLQi indicates correct normalization. In practice, the "correctly normalized" terms are 
normalized only at the tree level. Nevertheless, the double ratios suffer from uncertainties 
only of order ^^(A/mQ)^, even though the action is matched only at the l/mq level and the 
currents are matched only at the l/nicmb level. 

Once one is content to neglect corrections of order as^^/rnqY' , it is easy to obtain the 
continuum normalization of the lattice currents. By comparing the heavy-quark expansions 
for and to those for and A^, one sees that 

V^; = Py..V;t (2.12) 
A% = pA^.A^,b. (2.13) 

apart from discretization effects discussed above. The p factors are 

py., = CvJC'^l (2.14) 
p^c. = CaJC^. (2.15) 



and they are known at the one-loop level |TB 



The matrix elements are obtained from three-point correlation functions. For the zero- 
recoil B ^ D, B* ^ D* and B ^ D* transitions the three-point function are, respectively. 



C^^^(t/,t„ti) = Y.{^\ODix,tf)^,-f,^biy,t,)OU0,ti)\0), (2.16) 

x,y 

C^*^^* {tf, t„ t,) = E(0|(^D. (a;, t/)^c74^b(l/, QOi, (0, t,)\0), (2.17) 

x,y 

C^^^* (tf, t„ t,) = E(0|(^D. {x,tf)^,jjj,^fb{y, t,)oU0, Um, (2.18) 



x,y 
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where 0^(*) and O^ji*) are interpolating operators for the B^*^ and D^*^ mesons. In C^*^^* 
the spins of the vector mesons are parallel, and in C^~'^* the spin of the D* lies in the 
j direction. These correlation functions are calculated by a Monte Carlo method, as usual 
in lattice QCD. In the limit of large time separations, the correlation functions become 

V 2m£) a/ IniB 

C^'-^\tfM) = zie4/^i^^£2i^£)e---(*-*')e---(*/-*^) + ■ ■ ■ , (2.20) 
C^-^'{tf,t,,U) = 2^/24/2 i£l^£^25|^e-™-(*-*^)e-™-(*/-^^ + ■ ■ ■ , (2.21) 
where ttibm and m^(*) are the masses of the and D^*^ mesons. The normalization 



factors y Zfj(,) /2mfj(*) are conventional; they cancel when forming the double ratios, so 
we do not need them. The correlation functions defined in Eqs. ( p.l6| )-( p.l8| ) are the only 
objects needed from the Monte Carlo. In practice we hold tj = and t/ = T/2 fixed and 
vary ts over the range for which the lowest-lying states dominate the correlation functions, 
as is needed for Eqs. ( p.l9| )- (|2.21| ) to hold. (T = Nra is the temporal length of the lattice.) 



From the correlation functions we form the following double ratios 

o ^ g^-^(0,t,T/2) C^-^(0,t,T/2) 

C^-^(0,t,T/2) C^-^(0,t,T/2)' ^' ' 

P ^ C^*-^'(0,t,T/2) C^*-^'(0,t,T/2) 

c'^'-D*(o,t,T/2) C^'-^*(0,t,T/2)' ^ ^ 

C^-^' (0, t, T/2) C^-^* (0, t, T/2) 
^^i(^) - ^D^D- (0, t, T/2) C^-«* (0, t, T/2) • 

Apart from renormalization factors, these ratios correspond to the continuum ratios 7^+, 
TZi, and TZai- In the window of time separations t and T/2 — t for which the lowest-lying 
states dominate, all convention-dependent normalization factors cancel in the double ratios, 
and the ratios reduce to 



Pv-b\jR+ = V7^+ = /i+(l), (2.25) 
PvcbJR'i = Jn'i = hi{l), (2.26) 



Pa^^^Ra, = ^TZa, = hA,{l), (2.27) 

where p\ = pA^bpAt'^/ PA'^'^Pa^''- In particular, note that the axial current double ratio does 
not yield /iAi(l) directly, but instead /iai(1), defined in Eq. ( |1.10|) . Once we have com- 



puted the left-hand sides of Eqs. (|2.25|) - (|2.27|) for several combinations of the heavy quark 



masses, we can fit the mass dependence to the form predicted by the heavy-quark expansions, 
Eqs. dTTD-dni)- 

To summarize this section, let us review the steps needed to obtain the physical form 
factor JFb^£,. (1): 

1. compute the three-point correlation functions and thence the ratios Ri, Rai] 
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2. multiply -\/R+ and with py/r]y, and ^/Ra^ with Pa/va-, to obtain /i_|_(l)/?7y, 

/ii(l)/r7v, and /iai(1)/57a; 



3. fit 1 — /i/?7 [where h/rj is /i_|_(l)/?7v, /ii(l)/r7v, or /^^^(l)/?)^] to the heavy-quark mass 
dependence expected from Eqs. ( |l.ll|) -( p!.13|) ; 

4. use the resulting £y, and £p (and associated l/mg terms) to reconstitute hA^iX) = 
^B^D'iX) via (the l/mg version of) Eq. (pTSp. 

As discussed above, with the lattice action, currents, and normalization conditions cho- 
sen above, we obtain /iAi(l) with uncertainties of order as(A/2mc)^ and A^/(2mQ)^ from 
matching, although the fitting procedure also yields estimates of three of the four l/m^ 
terms in /iai(1), as discussed in Sec. |V|. 



III. LATTICE CALCULATION 

This work uses three ensembles of lattice gauge field configurations, which have been used 
in previous work on heavy-light decay constants [^,^], B ttIu and D — > ttIi^ semi-leptonic 
form factors ||2S|, light-quark masses [^, and quarkonia The quark propagators are 



the same as in Ref. ^7j, but we now use 200 instead of 100 configurations on the finest 
lattice (with (3 = 6.1). The input parameters for these fields are in Table |, together with 
some elementary output parameters. 



The quark propagators are computed from the Sheikholeslami-Wohlert (SW) action ^1 
which includes a dimension- five interaction with coupling csw, sometimes called the "clover" 
coupling. For the light spectator quark we use customary normalization conditions for 
massless quarks with the SW action, so csw is adjusted to reduce the leading lattice-spacing 
effect of Wilson fermions. In practice, we adjust csw to the value Mq ^ suggested by tadpole- 
improved, tree- level perturbation theory [^, and the so-called mean link uq is calculated 
from the plaquette. The leading light-quark cutoff effect is then of order agAa, multiplied 
by a numerical coefficient that is known to be small. For the heavy quarks we adjust csw to 
the same value, but, as explained in Sec. 0, one should think of this adjustment as tuning 
a coefficient in the HQET effective Lagrangian. 

The hopping parameter k is related to the bare quark mass. For the heavy quarks, Kh 
is varied over a wide range encompassing charm and bottom. For the light spectator quark, 
the first row of Hq in Table | corresponds to the strange quark. To test the dependence of 
the form factors on the light quark mass, we repeat the analysis for a few lighter spectator 
quarks. Table | also lists the tadpole-improved bare quark mass in GeV, 



amo = — — - , (3.1) 

Uo \2k 2/tcrit/ 

where the critical quark hopping parameter Kent niakes the pion massless. Although this 
mass is just a bare mass, it shows that the heavy quarks are heavy, and the hght quarks light. 

The lattice spacing a plays a minor role in our analysis, because both the lattice perturba- 
tion theory and the fitting to the heavy-quark mass dependence can be carried out in lattice 
units. The physical scale enters only in adjusting the heavy-quark hopping parameters to 
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TABLE I. Input parameters to the numerical lattice calculations, together with some elemen- 



tary output parameters. Error bars on the outputs 


TCiTCiT 'ir'i i" \r\t^ qcI" /Hi m "i" 1 c 1 
icicl LU Llic idoL U-i^lllo). 






Input; 






/? = 


6.1 


5.9 


5.7 


Volume, iV| x Nt 


X 48 


16^ X 32 


12^ X 24 


Configurations 


200 


350 


300 




1.46 


1.50 


1.57 


Kh, mo (GeV) 


0.080, 7.90 


0.077, 6.03 


0.062, 6.16 




0.090, 5.82 


0.088, 4.36 


0.089, 2.87 




0.097, 4.62 


0.099, 3.06 


0.100, 2.03 




0.100, 4.16 


0.110, 2.02 


0.110, 1.42 




0.115, 2.21 


0.121, 1.16 


0.119, 0.96 




0.122, 1.46 


0.126, 0.83 


0.125, 0.69 




0.125, 1.16 






Kg, rriQ (GeV) 


0.1373, 0.092 
0.1379, 0.039 


0.1385, 0.088 
0.1388, 0.073 
0.1391, 0.057 


0.1405, 0.09: 


t range 


[9, 15] 


[6, 10] 


[4, 8] 


Elementary outputs 


^crit 


0.13847t^ 


0.14017+^ 


0.14327i^ 


«rp-is (GeV) 
aj! (GeV) 


2-64+1^ 
2.40li° 


1.8lt^ 
1.47;^ 


0.89t^ 


Uo 


0.8816 


0.8734 


0.8608 


ay(3.40/a) 


0.14533 


0.15938 


0.18265 



the physical mass spectra, and in studying the dependence of /iAi(l) on a. Table | contains 
two estimates of the lattice spacing, from the spin-averaged IP-IS splitting of charmonium, 
Amip_is, and from the pion decay constant 



The renormalized strong coupling 0^/(3.40/0) at scale 3.40/a is determined as in Ref. p2 



In Sec. fV] the couphng is run to av{q*), where q* is the optimal scale according to the 



Brodsky-Lepage-Mackenzie (BLM) prescription p3|,p2|. Then av{q*) is used to calculate 



the short-distance coefficients pv/vv and Pa/va, which are introduced in Eqs. ( |2.25|) -( p.27|) , 
as well as the coefficient rjA- 

The right-hand side of Eq. (|2.19|) is the first term in a series, with additional terms for 
each radial excitation [and similarly for Eqs. ( |2.2CI| ) and ( p.21| )]. We reduce contamination 
from excited states in two ways. First, we keep the three points of the three-point function 
well separated in (Euclidean) time. The initial-state meson creation operator is always at 
ti = and the final-state meson annihilation operator at tj = Nt/2. We then vary the time 
ts of the current, to see when the lowest-lying states dominate. The second way to isolate 
the lowest-lying states is to choose creation operators and annihilation operators O^i*) 
to provide a large overlap with the desired state. This is done by smearing out the quark 
and anti-quark with IS and 2S Coulomb-gauge wave functions, as in Ref. [pl[] . 

Figure |l] shows the isolation of the ground state in the ratios R+{t), Ri{t), and RAi{t). 
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In each of the three modes we find a long plateau. We fit to a constant and obtain a 
precision at the percent level. For each ensemble, we choose the same fit range for all mass 
combinations listed in Table |. In Fig. |I] the resulting central values and error envelopes are 
given by the solid and dotted lines, respectively. Different fit ranges lead to slightly different, 
though consistent, results; this variation is folded in with the statistical error. Statistical 
errors, including the full correlation matrix in all fits, are determined from 1000 bootstrap 
samples for each ensemble. The bootstrap procedure is repeated with the same sequence 
for all quark mass combinations, and in this way the fully correlated statistical errors are 
propagated through all stages of the analysis. 

Figure ^ also demonstrates a clear distinction between the B* D* and the other two 
modes. Consequently, one can already see that iy is definitely greater than ip and £a, as 
expected from Refs. |TO|- p!^] . This is an important observation, because the largest l/rriQ 
correction to ^yii(l) is £y/(2mc)^. 



IV. PERTURBATION THEORY 



In this paper perturbation theory is needed to calculate the short- distance coefficients 
Pj (J = V, A) defined in Eqs. ( |2.14| ) and ( |2.15D , and rjj and 17^ appearing in Eqs. ( |1.5[ ) 
and ( |l.ll| )-( p!.13| ). The p factors match lattice gauge theory to QCD, and the t] factors 
match HQET to QCD. To fit the heavy-quark mass dependence of the lattice double ratios, 
one must also match lattice gauge theory to HQET, and the corresponding factors are 
simply pv/vv and pa/va- Figure | illustrates how these matching factors connect lattice 
gauge theory and HQET to QCD, and to each other. 

Lattice perturbation theory often yields a series that appears to converge slowly. The 
two main causes of the poor convergence have been identified the bare gauge coupling 
is an especially poor expansion parameter, and when tadpole diagrams occur expansion 
coefficients are large. These two problems can be avoided by using a renormalized coupling 
as the expansion parameter and by using perturbation theory only for quantities in which 
tadpole diagrams largely cancel. Then lattice perturbation theory seems to converge as well 
as perturbation theory in continuum QCD. 

To calculate the p factors only the vertex function is needed. By construction the self- 
energy contribution to wave-function renormalization, in particular the tadpole diagrams, 
cancels completely. Furthermore, even the vertex functions cancel partially, so the expansion 
coefficients should be small, as verified explicitly at the one-loop level [0. Indeed, as 
mqa —>■ 0, p —>■ 1, and as mqa —>■ 00, p —>■ rj. Thus, despite the fact that only the one-loop 
correction to pj is available 



it seems likely that perturbation theory can be expected 
to behave well, especially when measured against other uncertainties in this calculation. 

The other ingredient needed for an accurate perturbation series is a suitable renormalized 
coupling. We use the coupling ay defined through the (Fourier transform of) the heavy quark 
potential, as suggested in Ref. [^]. The scale q* of the running coupling av{q*) is chosen 



according to the BLM prescription [23,22 



log(g*a)^ 



(4.1) 



12 



1.04 



1.02 
1.00 
R 0.98 
0.96 
0.94 
0.92 



t 1 







♦ 

—1 

10 



12 



1.04 I : 1 : 1 : 1 : i : 1 : i : 1 r 




16 



(b) t 



1.04 I 1 \ 1 1 1 \ 1 \ 1 1 1 \ 1 \ 1 \ 1 \ 1 1 1 \ r 




FIG. 1. Double ratios RAi{t) (triangles), R+{t) (diamonds), and Ri{t) (squares) at (a) /3 = 5.7, 
(b) f3 = 5.9, (c) f3 = 6.1. The heavy quark hopping parameters are (a) {ki,,Kc) = (0.062,0.100), 
(b) {Kf„Kc) = (0.088,0.121), and (c) {Kb,Hc) = (0.097,0.122). The light quark mass is close to the 
strange quark mass. The lines represent constant fits in the indicated ranges. 
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ph 



QCD 



HQET 



FIG. 2. Diagram illustrating how the matching factors rj, and p/ij match lattice gauge theory 
and HQET to QCD, and to each other. 



where ( is pv/w or pa/va when fitting the mass dependence of the double ratios, or rjA 
when reconstituting /iAi(l) with Eq. ( p.. 51 ). The numerator *(^^'^ in Eq. ( ^.11 ) is obtained from 
the Feynman integrand for ^[^1 by replacing the gluon propagator D(k) by log{k'^a^)D{k), 
where k is the gluon's momentum. Such terms arise at the higher-loop level, so the BLM 
prescription sums a class of higher-order corrections. Since in the cases at hand the one- 
loop integrals are ultraviolet and infrared finite, the only scales that can appear are ^Jmjnh 
and 1/a. In general we find g* to be a few GeV; the only exceptions occur when (py/w)'^^ 
or [paI^a)^^^ are accidentally very small. 

One of the advantages of the BLM prescription is that the scale depends on the renor- 
malization scheme, in such a way that the value of the coupling itself does not depend on 
the scheme much. The coupling in an arbitrary scheme S is related to the V scheme by 



(47r)2 (47r)^ 



9W) gW) 



(0) 



(4.2) 



where for light quarks /5o = 11 



2nf/?), and 6^''' is independent of Uf. In many cases. 



the term dominates; for example, for the MS scheme, = —5/3 and 6^ - 

one chooses = g*e~*s'^/^, then g^igs) differs from Qyiq*) only by "non-BLM" terms of 
order 5'^(/9ofl'^)'~^) ^ ^ 2, which often are not very important. 
In summary, we evaluate all short- distance coefficients with 



l(0) 



-8. If 



[1] 



(4.3) 



and the appropriate BLM scale q*. To check for the possible size of non-BLM two- loop 
corrections (which are unavailable for pj), we also perform cross checks with aMs(9'lis)' 
obtain av{q*) via two-loop running from p2 



av'(3.40/a) 



2ai 



xl 



1 + VI -4.74aixi 



(4.4) 



where aixi = —{S/tt) Iiiuq. Uq and ay(3.40/a) are tabulated in Table |. 

Table |n| contains the values of pv/vv and Pa/va appropriate to the heavy quark mass 
combinations used in Sec. |V|. As expected, the perturbative corrections to these factors 
are small. The lattice coefficients and were obtained in Ref. 



18 . The continuum 



coefficients are p2| 
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TABLE II. Double ratios, computed in the Monte Carlo, and (re) normalization factors, com- 
puted in perturbation theory to one-loop BLM order. 



/3, Kq 










Pvhv 




PAlf]A 


6.1 


(0.080, 


0.115) 


1.00101^5 


1 7/1 

0.9851^^7 


1.0021 


1.0024+^^ 


0.9940 


0.1373 


(0.080, 
(0.090, 


0.122) 
0.100) 


1.0030l|[!^ 
l.OOOl^j]^ 


0.9742+1?^ 
0.9990+[]^ 


1.0008 
1.0002 


1.0043t?^^ 
1.00021°^ 


0.9919 
1.0000 




(0.090, 


0.125) 


1.0050^^? 


0.9757^1^ 


0.9978 


1.005ll^^ 


0.9908 




(0.097, 


0.115) 


1.0007^^^ 


0.9948l|} 


1.0003 


1.0012lJ^ 


0.9985 




(0.097, 


0.122) 


1.0023li 


0.987ll^^ 


0.9991 


1.0027l|| 


0.9954 




(0.100, 


0.125) 


1.0039+1^ 


0.9838+^^ 


0.9973 


1.00341^^ 


0.9933 


5.9 


(0.077, 


0.110) 


1 ^ A 

0.99811^1 


0.98721^;] 


1.0030 


1.0009t^^ 


1.0001 


0.1385 


(0.077, 
(0.077, 


0.121) 
0.126) 


0.9971+^^ 
0.9984+^^ 


0.9697+^J 

1 fin 

0.95491^? 


1.0035 
1.0015 


1.0030l^^ 

1 '7n 

1.00541^0 


0.9770 
0.9868 




(0.088, 


0.110) 


0.9993ljf 


0.9934lj| 


1.0013 


1.0007lJ^ 


0.9999 




(0.088, 
(0.088, 


0.121) 
0.126) 


0.9993+^^ 
l.OOlll^^ 


0.9795+1^ 
0.9666^^? 


1.0016 
0.9995 


1.00281P 
1.0053138 


0.9944 
0.9903 




(0.099, 


0.110) 


0.9999l[j^ 


0.9980l[j^ 


1.0003 


1.0003l|]^ 


0.9990 




(0.099, 


0.121) 


1.0003^1^ 


0.9883^1^ 


1.0000 


1.0019lJ^ 


0.9969 




(0.099, 


0.126) 


1.0022^11 


0.9780l^^ 


0.9983 


1.004ll|^ 


0.9983 


5.7 


(0.062, 


0.089) 


Tin 

0.99441^^ 


0.99231^^ 


1.0024 


0.99751^^ 


1.0010 


0.1405 


(0.062, 
(V\ 069 


0.100) 
n 1 95^1 


0.9895+^^ 
9786+^''^ 


0.9845+^^ 


1.0050 
1 01 1 4 


0.9958t^^ 

n 9888+-'^^-'^ 


1.0017 




(0.089, 


0.100) 


0.9992+[]^ 


0.9984l[]^ 


1.0005 


0.9996t|]^ 


1.0001 




(0.089, 


0.110) 


0.9969^^^ 


0.9929^11 


1.0018 


0.9985lJi 


1.0002 




(0.089, 


0.119) 


0.9945li 


0.9816li 


1.0035 


0.9969111 


1.0000 




(0.089, 


0.125) 


0.9939l|^ 


0.9673^^^ 


1.0041 


0.9958l|^ 


1.0112 




(0.100, 


0.125) 


0.9979lJ^ 


0.9793^^^ 


1.0022 


0.99831^? 


0.9958 




(0.110, 


0.119) 


0.9997^[j^ 


0.9972l[j^ 


1.0004 


0.9998l[J^ 


0.9995 



r^W = Cp3/(m5/m,)/167r2, (4.5) 

1]^ = Cp 9/(mb/mc)/327r^ + ijy^ \n{mbamca), (4.6) 

r/S^ =C^3/(m,/m,)/167r2, (4.7) 

Cp 5/(m5/mc)/327r^ + 77^4' \n{mbamca) , (4.8) 



Va 



where 



/(2) = i±ilnz-2. (4.9) 



The important properties of f{z) are /(I) = 0, f{l/z) = f{z). From the matching procedure 
derived in Ref. [jl8[ one sees that the masses used in f{mi,/mc) should be the kinetic masses. 



namely the mass appearing in the kinetic term in Eq. ( |2.1| ). 

Two different schemes for defining the kinetic quark mass are used in this paper, because 



they are simple to implement. Both employ the formula |15 
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+ (4.10) 



am2 C^*"! sinh(ami) e 

which is the tree-level relation between the kinetic mass m2 and the rest mass mi, for the 
SW action. One choice is to use the tree-level value for the rest mass ami = log(l + amo), 
with amo from Eq. ( |3.1| ), and we call the result the tree- level kinetic mass. The other is 
to use the one- loop rest mass in Eq. ( |4.10|) p3[| , and we call the result the quasi-one- loop 
kinetic mass. (The kinetic mass receives further radiative corrections, but they are known 
to be small [^].) The second choice is essentially the (one-loop) perturbative pole mass. 
Although the difference between these schemes is formally of the non-BLM two-loop order, 
they could give slightly different results in practice. Thus, using both and comparing gives 
us a handle on the terms omitted from the perturbative series. 

When reconstituting the physical form factor /i^i (1) with Eq. ( |1.5| ), one needs a numerical 
value for the short- distance coefficient rjA- Although it is known at the two-loop level 
we use the one-loop, BLM results, so that all perturbation theory is treated on the same 
footing. Thus, we take |^ 



r/W = Cf [3/(m5/m,) - 2] /IGti', (4.11) 



*Va = Cf [lf{mb/mc) - ij /IQn' + r]'^^ \n{mham,a). (4.12) 

For consistency, it is necessary to use the same definition of the quark mass in rjA as in p/r]. 

If we take the quasi-one-loop kinetic masses, which are very close to continuum pole 
masses, we find z = m^dm^y, = {0.308,0.296,0.290}, q* = {2.94,3.08,3.12} GeV, aviq*) = 
{0.205,0.203,0.208} and, hence, 

r/^ = {0.9713, 0.9724, 0.9724} (4.13) 

for (3 = {5.7,5.9,6.1}, respectively. On the other hand, if we take the tree-level ki- 
netic masses, we find z = {0.221,0.230,0.234}, q* = {2.02,2.14,2.14} GeV, ayiq*) = 
{0.241,0.238,0.245} and, hence, 

r/^ = {0.9769, 0.9758, 0.9746} (4.14) 

for j3 = {5.7, 5.9, 6.1}, respectively. Note that although the coupling is larger in this scheme 
(because the quark masses and, hence, q* are smaller), the perturbative correction is smaller, 
because the magnitude of the coefficient 77^4' decreases with z. As we shall see below, this 
scheme dependence in rjA is largely cancelled by the corresponding scheme dependence of 
the l/mg corrections. 

These values of rjA are slightly larger than the value 0.960 [^,0, which is widely adopted 
in the literature. The origin of this difference is the value used for a^. We extract from 
lattice QCD, which, in the quenched approximation, underestimates slightly ||3^. Also, 



[1] 



there is nothing special about the standard value. It does not include uncertainties from the 
measured value of as{Mz) or from the b and c masses. When our method is applied to full 
QCD, the double ratios, the gauge coupling, and the quark masses all can be determined 
self-consistently. In the meantime, we shall assign uncertainties from omitting the non-BLM 
two-loop term, adjusting the heavy quark masses, and the quenching effect on a^. 
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V. HEAVY QUARK MASS DEPENDENCE 



In this section we fit the (suitably normahzed) double ratios to the form expected from 
the heavy quark expansion, yielding the quantities a^iy, cl^^a, and a^ip (i.e., in lattice 
units). We find that it is also necessary and beneficial to incorporate terms of order l/rriQ 
in the heavy quark expansion. The last step is then to combine these results into the main 
goal, which is hA^il)- 

Table |T| contains the results of our Monte Carlo calculations of 



addition to the short-distance coefficients discussed in Sec. This information is combined 
to form 

(5.1) 



Vv Vv 
Pv^/R^ _ hi_ 

Vv Vv' 

PaJRai Kai 



(5.2) 

- , (5.3) 
Va Va 

which we fit to the expected heavy-quark mass dependence. For each ratio in Eqs. (|5.1|) - (|5.3|) 
we try the fit 

^ = l-iA^c(^) + lc(3)S,), (5.4) 
where S-"^^ and S-^^ are taken as free fit parameters, and 

^2 = — (5.5) 

am2c CiTn2b 

S2 = + . (5.6) 

am2c ci'm2b 

In Eqs. (|5.5| ) and (|5.6| ), the subscript 2 indicates that the kinetic mass m2 appears. For 
the quadratic term we use Ag, even though the masses m2, mg, and all appear in the 
heavy-quark expansion to lattice QCD EO] , because m2 = mis = f^s at our level of accuracy. 



The rest mass mi in Eq. (2J.) drops out completely |20 . 

The l/rriQ term is introduced in Eq. ( |5.4| ) to describe the data over a wide range of I/uiq. 
The particular form A^S is the only one that is invariant under the interchange symmetry 
c b and vanishes for rric = mb- The l/rriQ terms arise from many sources in HQET. Some 
of them, like triple insertions of the I/ttlq terms in >Chqet, are correctly normalized with 
the choice of lattice action and currents made in Sec. |T|. They lead to A2S2, with (to our 
accuracy) the kinetic mass everywhere. Others, like an insertion of a l/rriQ term combined 
with an insertion of a l/mq term, are not and would lead to A2Ax^x, where Ax'^x amounts 
to the difference of short- distance coefficients for the higher-dimension HQET operator Ox- 

The most important mismatches of Ax'^x are of order asam2c and of order (am2c)^, 
provided am2c < 1. They are not necessarily small but, perhaps, small enough to pin down 
the l/rriQ corrections. The l/rriQ contributions are influenced mostly by the region with large 
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S, where am2c < 0.6. Thus, the fit coefficients c^^^ can be expected to give a reasonable 
estimate of the desired a^i^^\ Moreover, corrections of order (K/mqY are small to begin 
with, so even a large relative uncertainty in them leads to a small absolute uncertainty 
on /lAi(l)- 

As mentioned in the Introduction, there are four l/mq terms in the heavy quark expan- 
sion of If we write 



1 + <^l/m2 + 



then can be read off by comparing with Eq. (|L 



and 



^(3) 



+ 



?(3) 



(2m,)3 (2m,)2(2mb) (2me)(2mb)2 (2mfe) 



(5.7) 



(5.8) 



As suggested by the notation, iy'^ is related to hi{l), and ip' is related to /i+(l). Repeating 
the argument based on heavy-quark spin symmetry, first for the b, then for the c, one sees 
that /tAi(l) and hi{l) share the term iy^ / {2mc)^ , and that /iAi(l) and /i+(l) share the term 
£?V(2m,)^ as given in Eq. — 



«(3) 



^8]). The other two terms in Si/m^ can be rewritten 



^(3) 



+ 



?(3) 



-A 



(2mcYi2mh) (2mc)(2m;,)2 (2mc)(2m;,) \2mc 2m;,/ (2mc)(2mb) \2m, 



+ 



+ 



(<{3) 



2mh 



(5.9) 



where 



^(3) 



7(3) 



/2 and ^^^^ 



D 



-B 



/2. Simple algebra shows that £^^^ is 
indeed the coefficient of the A^S term in the heavy-quark expansion of the ratio /^^^(l). 
Thus, to the extent that we can identify C|py^| with a^^jpy^}, we can reconstruct three of 

the four l/rriQ corrections to hAi{l)- Only £^^^ eludes us. 

To show the quality of the fit to the mass dependence, we plot in Fig. ^ the quantity 



Q 



1 - pVR/t] 



lc(2) + lc(3)S, 



(5.10) 



vs. S2 = (l/am2c + l/am2b), with the quasi-one-loop definition of am2. Linear behavior 
in (l/am2c + l/a'm2b) is observed for each form factor, and we show the fit line in the 
figure. Some curvature is noticeable for the heaviest masses in Fig. ^(a), but the linear 
fit is still consistent within statistical errors. The growth of the statistical error toward 
the heavy-quark limit is a property of the heavy-light meson in the Monte Carlo, and it is 
unavoidable [l35|,B6i. 



The values of the fit parameters cjp 



V,A} 



(3) IT -n 

{p,v,A} ^^"i '^{PVA} listed in Table |T] 



In each case the extracted values of c^^-* and c^^^ are highly correlated. On the other hand, 
the combinations 



,.(2) I 1 J3 
<^A + o^A 



S3) 



2am2c 
(3) f 



+ 



\2amor 2am 



2b 



.(2) 



+ 



2 am 



2b 



(5.11) 
(5.12) 

(5.13) 



18 



0.32 




-0.08 



1/am, + l/am^. 

2c 2b 



0.08 




(b) 



1/am, + 1/am,, 

2c 2b 




(c) 



+ l/flm„, 

2c 2b 



FIG. 3. (1 — h/r])/A2 vs. l/am2c + l/am2b when /i/r/ is /ii(l)/r?y (squares), /i_|_(l)/r/v' (dia- 
monds), and hA^{l)/fiA (triangles) at (a) (3 = 5.7, (b) /? = 5.9, (c) (3 = 6.1. Solid lines are best fits 
and dotted lines are error envelopes. 
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TABLE III. Coefficients in tlie l/mg expansion, Eq. (^). 



p 


h, 

Cp 


-/w 

Cp 




(2) 


4^ 


J2) 


hAi/fjA 


J3) 
^A 


6.1 


-0.0191°^^ 


0.015 


f035 
-035 


0.117;°! - 




0.0351 


050 
042 


0.018l[]i 


5.9 


-0.0141°^ 


0.012" 


f033 
-030 


0.184;°^^ - 


0.089;oi 


0.0371 


040 
042 


0.022l°^2 


5.7 


0.075t?^° 


-0.100" 


f099 
-123 


0.289;j^^ 


0.145;ii 


0.0891 


099 
118 


0.03011^^ 



are statistically more precise, because the correlated error cancels, for the first two especially 
so. These combinations appear directly in Eq. ( |5.7| ), provided we can reliably identify Cy 
with a^Vy . We argued above that this identification is not too bad, because the coefficients 
c^^-* should be influenced principally by smaller masses. As seen in Fig. ^, this predjudice is 
borne out, especially when the correlated statistics are taken into account: the best fits fit 
best for large {l/am2c + l/am2b)- 

The results presented in Fig. ^ and Table |T| are all for the quasi-one-loop definition 
of am2- One should keep in mind that the is and have a well-defined interpretation 
as matrix elements within HQET. Their detailed definition depends on the renormalization 



scheme of operators in HQET, as discussed, for example, in Ref. |^^. After reconstituting 
/lAi(l), however, the scheme chosen should have only a minor, residual effect. Repeating the 
fits with the tree-level definition of changes the fit coefficients significantly (as expected). 
The change in /iAi(l) is, however, not great, and it is of order a^/mg, as expected. 

To fix the physical values of mi, and we compute the Bg and Dg spectra on the same 



ensembles of lattice gauge fields. Combining these inputs with the second row of Table |T] 
{(3 = 5.9) (and omitting the contribution) we find 

^l/m" = 5l/m2 + Si/jn3 (5-14) 

fed o/JefT feS 

I _ n Q 4 47+0.0078 \ /c-.c-N 



(2mc)2 2mc2mb (2m 



b 



which is needed in Eq. (|5.7|) . The error quoted here is statistical only; systematic uncer- 
tainties are considered in detail in the next section. Equation ( |5.15| ) shows the power of our 
method: even with 15% statistical uncertainties on 6i/m" = hAi/VA — 1, one can see that 
/i^j(l) itself can be very precise. 



VI. SYSTEMATIC ERRORS 

The intermediate result in Eq. (|5.15| ) is obtained at one value of the lattice spacing, and 
with a spectator quark whose mass is close to that of the strange quark. In this section 
we consider the systematic uncertainty from varying a and rriq, as well as those from other 
sources. Table |V| summarizes the results of this analysis, giving the absolute error on the 
main result, /i^j(l), and also fractional error on 1 — /i^j(l). As noted above, the uncertainties 
should scale with 1 — /i^j(l). 

In the following subsections, we consider, in turn, the uncertainties arising from fitting 
Ansatze, which incorporate contamination in Eqs. (|2.19| )-( p.21| ) of excited states (Sec. |V1 A| ); 
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TABLE IV. Budget of statistical and systematic uncertainties for hAi{l) and 1 — The 
row labeled "total systematic" does not include uncertainty from fitting, which is lumped with the 
statistical error. The statistical error is that after chiral extrapolation. 



uncertainty 










1 - /iAi 
(%) 




statistics and fitting 


+0.0238 




-0.0173 


+27 




-20 


adjusting ttIc and 


+0.0066 




-0.0068 


+ 8 




- 8 






+0 0082 






± 9 




a,(A/2mQ)2 
(A)V(2mQ)3 




±0.0114 






±13 






±0.0017 






± 2 




a dependence 


+0.0032 




-0.0141 


+ 4 




-16 


chiral 


+0.0000 




-0.0163 


+ 




-19 


quenching 


+0.0061 




-0.0143 


+ 7 




-16 


total systematic 


+0.0171 




-0.0302 


+20 




-35 


total (stat © syst) 


+0.0293 




-0.0349 


+34 




-40 



heavy quark mass dependence (Sec. VI B| ); matching lattice gauge theory to HQET and 
QCD (Sec. |VI C|) ; lattice spacing dependence (Sec. [VI D|) ; light (spectator) quark mass 
dependence (Sec. |VI E] ); and the quenched approximation (Sec. [VI F]) . In Table |I^ the 
statistical uncertainty is added in quadrature to that from fitting, as discussed in Sec. [VIA. 



As outlined in Sec. |T|, statistical uncertainties are computed with the bootstrap method 
and full covariance matrices. 



A. Fitting and excited states 

We define iii our fits with the full covariance matrix. For the plateau fits to R{t) 

= E [^(^i) - ^fit] cr-\tut2) [Rit2) - Rfit] ■ (6.1) 



tl,t2 



Because the numerical data are so highly correlated, some components of the (inverse) 
matrix cr~^(ti,t2) cannot be determined well. These components are discarded, according 
to singular value decomposition (SVD), by eliminating eigenvectors of o"^ whose eigenvalue 
A < TsvoAmax, with rgyo Small. We find we have to set rsvo ~ 10~^ to remove the noisy 
eigenvectors from in Eq. ( [b.l[ ). 

A potential drawback of the double ratio technique is that an early plateau could be 
induced. We cope with this issue by trying many fit ranges for the time tg of the current. In 
general, fits to a constant have good ^-^^d agree for fit ranges within the plateaus clearly 
seen in Fig. |I]. For each ensemble of lattice gauge fields we choose a single range for tg for 
all three ratios and all heavy quark mass combinations. In each case, the range is chosen to 
give small statistical error on Rat, while maintaining a central value close to that from short 
intervals centered on T/4. 

The expressions in Eqs. ( p.l9[) -( ^.21[ ), relating three-point correlation functions to matrix 



elements, suppress terms from radial excitations of the desired, lowest-lying states. Because 
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of heavy-quark symmetry, corresponding excitations of the D and B systems have similar 
wave functions and mass sphttings. Consequently, their contribution to the double ratios 
largely cancels, leaving a residue that is suppressed by (A/mg)^ as well as the exponential 
factor for large times. Thus, the excited-state contamination in a double ratio scales as 
R — 1, rather than R. 

The fits of the heavy quark mass dependence are obtained by minimizing 

= E {Q^ - jc^'^ - -sc^''^2.) arf (Q, - Ic^^) - Ic^^)^^,) , (6.2) 

where i, j label mass combinations. Once again, not all components of a'"^ are well deter- 
mined. The fits are stable with rgyo = {5 x IQ-^, 5 x 10-^ 1 x IQ-^} for (3 = {5.7, 5.9, 6.1}. 

In summary, the fitting procedure to determine the double ratios R+, Ri, and Rai 
depends on the fit range for ts and on the cut rsvo in the SVD. Similarly, the fit parameters 
of the heavy quark mass dependence, c^^-* and c^^\ depend on an additional SVD cut. The 
central values quoted here are from the fit ranges given in Table |, rsvo = 10~^ for R{t), and 
rsvD as given above for c*^^-* and c^^\ We then repeat the analysis with larger and smaller 
SVD cuts and, for R{t), with other fit ranges. The resulting variation in /iAi(l) is smaller 
than the statistical error of the "best fits" . Since excited states contribute differently in each 
fitting Ansatz, the uncertainty in fitting R{t) incorporates the uncertainty due to excited- 
state contamination. For convenience in analyzing the other systematics, the fitting error is 
added in quadrature to the statistical error. 



B. Heavy quark mass dependence 

The physical heavy quark masses enter when reconstituting with Eq. (|5.7| ). We 
determine them by tuning the hopping parameters Kb and Kc to reproduce the Bs and Dg 
spectra. To do so, we must compute the meson kinetic masses, which are somewhat noisy, 
and we must choose an observable to define the (inverse) lattice spacing. Thus, the tuned 
values of Hf, and Kc have statistical uncertainties, from both the meson masses and from a~^. 

They also have systematic uncertainties. For example, the inverse lattice spacing 
is not the same when defined by the IP- IS splitting of charmonium or by /t^, as noted in 
Table |. Similarly, Kb and Kc are not the same when quarkonium spectra are used instead 
of heavy-light spectra, although for Kc this makes very little difference. In the end, we are 
left with a range for Kb and k^ and, hence, the heavy quark masses used in Eq. (|5.7|) . This 
range leads to the error bar labeled "adjusting and m^" in Table [TV] . 



C. Matching 

As discussed in Sec. || our method for heavy quarks matches lattice gauge theory to QCD 
by normalizing the first few terms in the heavy-quark expansion This is necessary 

to keep heavy-quark discretization effects under control, but the approximate nature of 
the (perturbative) matching calculations leads to a series of uncertainties. The three most 
important of these are listed in Table fV. 
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TABLE V. Scheme dependence of /lyi^(l). For each value of /5, the columns compare the 
scheme with tree-level and quasi-one-loop kinetic masses in r]j and in the mass dependence. The 
rows compare the effect of the contributions; here S = l/rric + l/mf, refers to the correction 



m Eq. Each row includes the corrections from all preceding rows. 



1/m" 


f3 = 


6.1 


P = 


5.9 


p = 


5.7 




tree 


quasi 


tree 


quasi 


tree 


quasi 


IMq 

T./{mcmb) 


U.OrOO_o.o372 

n Qqqn+0.0150 
U.yOOZ_Q Qj^24 

u.yz/o_o 0114 


0.93291°:°?°? 
0.93261°:°?°^ 
0.9274l°:°}i 


0.9078;°:°^^? 

0.93621°:°°^? 
0.9363l°:°°i? 

n Qqqo+0.0057 
U-y<3 -30-0.0052 


0.9103l°:°J^° 

n QQOi +0.0082 
u.yo/:i_Q QQ72 

n qq9n+0-0082 
qS00+°-°°^^ 


0.9365l°:°ll? 
0.9549l°:r 

n qc^fil +0-0099 
U.yoOi_o.o086 

U.yOUO_o.0079 


0.9303l°:°237^ 
q4n6+° °i^2 

U.y4UD_Q Q^5^ 

0.9409l°:°163 
q4n0+°-°^^2 

U.y4UU_Q Q135 



The first is formally of order a^. It comes from omitting the non-BLM radiative correc- 
tions to the factors pj and t]j and from omitted loop corrections to the quark masses and 
to ag. As discussed in Sec. 0, pj comes from the cancellation of (continuum and lattice) 
vertex functions. Thus, by design, the coefficients of its perturbation series are small — 
usually smaller than those in rjj |T3. With rjA (and rjv) we can check explicitly how big the 



non-BLM two-loop corrections are. For example, the value of /iAi(l) is reduced by 0.0082 
if we switch to the MS scheme and include the non-BLM two-loop part of the r]j. Since 
the unknown two-loop corrections to the pj could compensate, or even over-compensate, we 
take the two-loop uncertainty to be ±0.0082. 

The next matching uncertainty is formally of order as{A/2mc)'^, from tuning the lattice 
action and currents to HQET. Setting = 0.2, A = 500 MeV, and = 1.25 GeV, one 
finds a,(A/2mc)2 = 0.008. Another way to estimate this effect is to compare the analysis 
with tree-level heavy quark masses to the standard one with quasi-one-loop masses. The 
difference in /2-Ai(1) is in the same ballpark, at most +0.0114. Since other schemes for the 
quark mass could lead to shifts in the other direction, we take ±0.0114 as the uncertainty 
from this source. 

The last matching uncertainty is of order (A/m^)^, from the omission of 



1 1 



0.0017, (6.3) 



(2mc)(2mfc) \2mc 2m{, 

assuming f^"* = A^, = 4 GeV, and the same values as above. With same choices made 
above, we estimate that / {2mc) + l/{2mb)]/{2mc2mb) and ip^ /{2mb)^ should be around 
0.0033, and 0.0002, respectively. In Table we show the actual effect of the included l/rriQ 
corrections. The scatter of the different analyses bears out the latter estimates, lending 
credence to Eq. (|6.3|) . Uncertainties in the included l/mg terms are smaller than Eq. ( |6.3|) , 
because many of them are obtained correctly, and the mismatch in the others is small. 



D. Lattice spacing dependence 

The lattice calculation of /i^i has lattice artifacts from heavy quarks, light quarks, and 
gluons. For the heavy quarks, discretization effects and heavy-quark effects are inevitably 
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FIG. 4. Lattice spacing dependence of /iAi(l) (triangles), (diamonds), and /ii(l) 

(squares). The light quark mass is close to the strange quark mass. The solid (dotted) lines 
represent best fits (error envelopes). 



intertwined [0,0 , and are mostly part and parcel of the matching uncertainties considered 
above. The light quarks suffer from discretization effects of order a^Aa and (Aa)^; the gluons 
of order (Aa)^. That being said, we can test for the magnitude of discretization effects, by 
comparing the analysis of Sec. for three lattice spacings. The results are plotted against a 
in Fig. H, which also contains results for /i+(l) and hi{l). The last two are much closer 
to 1 and their statistical uncertainties are correspondingly smaller. This underscores, once 
again, that the uncertainties scale as 1 — h. 

The result for /iAi(l) with the available I/itlq contributions (solid triangles) is consistent 
with a constant. We take as our central value the average from the two finer lattices, because 
for them the (heavy-quark) discretization effects are smaller. This is 

hAAl)=0m93t'Zl2 (6.4) 

where the error is the statistical error on the average, with the error from fitting added in 
quadrature. In Fig. ^ the solid and dotted lines indicate this average and error band. 

The third point, at a = 0.84 GeV (from f3 = 5.7), has the greatest uncertainty from 
heavy quark discretization effects, so it is excluded from the central value. Instead we use it 
to estimate discretization uncertainties. If one assumes that discretization effects from the 
light spectator quark and gluons are negligible, then it would be appropriate to average all 
three. This average is slightly higher, and we take this increase as the upward systematic 
error bar. If, on the other hand, one assumes that the light spectator quark's discretization 
effects are responsible for the somewhat larger value of /iAi(l) on the coarsest lattice, then 
it would be appropriate to extrapolate linearly in a. The dashed line in Fig. ^ shows this 
extrapolation. The extrapolated value is significantly lower, and we take this decrease as 
the downward systematic error bar. The error bar resulting from these two estimates is very 
asymmetric: laoili- 
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1.00 



0.96 




0.04 0.08 2 0.12 0.16 

(m a) 

FIG. 5. Dependence of /lAi(l) at /3 = 5.9 on the mass of the hght spectator quark. Here 
is the mass of the pseudoscalar meson consisting of two "hght" quarks. The sohd (dotted) hnes 
represent the best hnear fit (error envelope). The lower (upper) curves with a cusp add to the 
linear behavior the contribution in Eq. (|6^), taking go'Dw = 0.60 {go'Dw = 0.27). 



E. Chiral extrapolation 



The calculations discussed so far have a spectator quark whose mass is near that of 
the strange quark. Figure || shows how /iAi(l) changes for lighter spectator quarks, on the 
lattice with (3 = 5.9, for which we have three values of the light quark mass, /iyii(l) is 
plotted against (in lattice units), which is a physical measure of the light quark mass. 
Since the statistical errors in Fig. |^ are highly correlated, the downward trend in /iAi(l) is 
significant. The same trend is seen for f3 = 6.1. Extrapolating linearly in to the physical 
pion, reduces the result in Eq. (|6.4|) to 



n qnn+0 0238 



(6.5) 



and increases the statistical error. This value, using the average of the /3 = 5.9 and 6.1 
lattices and the chiral extrapolation from /3 = 5.9, gives the central value in Eq. ( 1.14 ). 

In the chiral expansion, the terms responsible for the linear behavior are formally of 
order A^m^/(2mc47r/^)^. Terms of order A^/(2mc 47r/^)^ are larger for the physical pion 
mass, but are comparable for our artificially large pion masses. Randall and Wise ||^ have 
computed the m° effect at one loop in chiral perturbation theory. They find 



{2mc 



(2mc 



+ 



9d*Dtt 



' AW ^ 



2 r 



In^ 



+ /(- 



(6.6) 



where 

is the D*-D-7T coupling, A^'^) 



2mj^ is the mass of the pseudoscalar meson with two strange quarks, gD*Dir 



mo* — ran = 142 MeV is the D*-D mass splitting, and 



aw/ 



TC. 



rjs). For go'D-K we consider the range 0.27-0.60, which encompasses 



estimates based on fits to experimental data {qd* d-k 



n 07+0.O6 
•J-^'-o.os 



quark models {go* d-k 
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0.38 0), quenched lattice QCD (c/d-Dtt = 0.30 ±0.16 H or gB-'Sn = 0.42 ±0.09 |T1), and 



the recent measurement of the D* width {qo'D-k = 0.59 ± 0.07 |42 



The chiral loop function f{x) has rather different behavior, depending on x. At x = — 1, 
which turns out to be the physical region (m^r ^ A there is a cusp, and the value of / 
becomes large: /(—I) ~ 11 whereas /(— x^J = /(— 0.2) = 0.53. To illustrate this behavior, 
we have shown in Fig. ^ the sum of the second term in Eq. ( |6.6|) and the linear fit. In the 
region where we have data, the term from Eq. ( |6.6| ) hardly varies, but near the physical limit, 
it bends the curve down. With the quoted range for gD*Dn, the decrease in /iyii(l) amounts 
to 0.0033-0.0163, coming mostly in the region where ~ A'^'^^ as shown in Fig. ^. In an 
unquenched calculation, one would add this contribution to /i^j(l). Because Qwd-k remains 
uncertain and because we are using the quenched approximation, we take it as an additional 
systematic uncertainty of t^oiea- This effect and the amplification of the statistical error 
together make the chiral extrapolation the largest source of uncertainty. 



F. Quenching 

An important limitation of our numerical value for /i^i(l) is that the gauge fields were 
generated in the quenched approximation. The quenched approximation omits the back- 
reaction of light quark loops on the gluons, and compensates the omission with a shift in 
the bare couplings. Two obvious consequences of quenching are that the coupling as runs 
incorrectly, and that pion loops [as in Eq. ( |6.6| )] are not correctly generated. 



Let us consider first the effect on the running coupling. The values for 7]a in Sec. |^ are 
obtained with the quenched coupling. If is corrected for quenching, it is larger and 
the short-distance coefficients are changed by —0.0050 for rjA and ±0.0032 for rjy. These 
changes both reduce /lyi^(l). 

For the pion-loop contribution we can look to comparisons of quenched and unquenched 
calculations of other matrix elements. Studies of the decay constants fs and fo show 



discrepancies on the order of 10% between quenched and (partly) unquenched QCD ^3] , ^ 
A form factor, which is the overlap of two wave functions, is presumably less sensitive to 
quenching than a decay constant, which is a wave function at the origin. So, one should not 
expect the quenching error here to be more than 10%. Even in the quenched approximation 
all three double ratios tend to unity in the heavy-quark symmetry limit. Thus, the quenching 
error, like all others, scales with i? — 1, rather than R. We therefore apply the estimate of 
10% to the long-distance part, ^i/m^, to obtain an error bar of ±0.0061. 

We estimate the total quenching uncertainty to be the sum of these two effects, or lo;o?43. 

G. Summary 

Combining Eq. (|6.5| ) with the error budget in Table fV|, we obtain 

h n — n Q1 'in+O-O^SS+0. 0156+0.0032+0. 0000+0.0061 (a7\ 

"'Al\^) — U.J10U_Q Q]^73_Q Q]^57_Q Q;^4;^_g Q]^g3„Q Q]^43, I " • ' J 

where the error bars are from statistics and fitting, adjusting the heavy quark masses and 
matching, lattice spacing dependence, light quark mass dependence, and the quenched ap- 
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proximation. (The uncertainties on the second through fifth rows of Table |I^ are added in 
quadrature.) Adding all systematics in quadrature, we obtain 



/^A,(i) = :fb^d^{i) = o.9i3ol°:°?ll°:°J^^. (6.8) 

Although we have considered all sources of systematic uncertainty, it is not possible to 
disentangle them completely. For example, the lattice spacing dependence is not completely 
separated from the HQET matching uncertainties, and the quenched approximation affects 
the chiral behavior, the adjustment of rric and mb, and, through as, the matching coefficients. 



VII. COMPARISON WITH OTHER METHODS 

In this section we compare our method, based on lattice gauge theory, with others existing 
in the literature. To do so, it is convenient to refer to Eq. ( p..5|) and discuss how the short- 
and long-distance contributions are evaluated. 

One approach, sometimes advertised as "model-independent" , is to estimate the h with 
the non-relativistic quark model [p|,p!0|. The more recent estimate [1^ takes 5i/m2 to be 



—0.055 ± 0.025 by covering a range of "all reasonable choices". Combining it with the 
two- loop calculation of r;^, one obtains 

J^B^D*{1) = 0.907 ± 0.007 ± 0.025 ± 0.017, (7.1) 

where the quoted uncertainties |1T0|J6[] are from perturbation theory, errors in the quark model 



estimate of the l/wig terms, and the omission of l/rng terms. Uncertainties from and 
the quark masses are not included. A fair criticism of this approach is that it does not pay 
close attention to scheme dependence of the long- and short-distance contributions. The 
standard (/i-independent) result for tja corresponds to renormalizing the operator insertions 
of HQET in the MS scheme. The quark model estimates, on the other hand, are presumably 
in some other scheme, so there is a possibility to over- or undercount the contribution at 
the interface of long and short distances. 

Another approach is based on a zero- recoil sum rule [p!T],p|. These authors prefer to 



introduce a concrete separation scale /i. In this scheme rj^ and the is depend explicitly 
on fi. The /i-dependent two-loop part of rjA is known A recent estimate of the zero- 
recoil form factor is |46[ 



J^B^D'il) = 0.89 ± 0.015 ± 0.025 ± 0.015 ± 0.025, (7.2) 

where the quoted uncertainties are from the unknown value of the kinetic energy /i^(/i), 
higher excitations with D* quantum numbers and energy E < mD* + fi, perturbation the- 
ory, and the omission of I/ttlq terms. We note that both /x^ and the excitation contribu- 
tion should, in this scheme, cancel the (/i/mg)^ part of f]A{fJ')- Since there is no model- 
independent method to calculate the excitation contribution (except unquenched lattice 
QCD), it is not clear how to implement this cancellation. 

As shown in Fig. |], our result Eq. ( |1.14| ) agrees with the previous results, within errors, 
and the quoted errors are of comparable size. Our result includes an estimate of three of 
the four l/rriq contributions. All three are subject to a QED correction of +0.007 [0. An 
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important feature of our method is that, even in the quenched approximation, we are able to 
separate long- and short-distance contributions self-consistently. Indeed, we have repeated 
the calculation with two different schemes for the heavy quark masses, and the results are 
the same. Furthermore, it is clear that moving terms of order /i^/mg between the long- and 
short- distance parts will cancel out in our method, as long as it is done consistently. Finally, 
with future unquenched calculations in lattice QCD, our method allows for a systematic 
reduction in the theoretical error on \Vcb\- 



VIII. CONCLUSIONS 



We have developed a method to calculate the zero recoil form factor of B —>■ D*!^ 
decay. We introduce three double ratios in which the bulk of statistical and systematic 
errors cancels, thus enabling a precise calculation of J-'b-,d*{^)- By matching lattice gauge 
theory to HQET, we are able to separate long-distance from short-distance contributions. 
Then the coefficients in the l/mq expansion are obtained by fitting the numerical data. In 
this way we obtain the (leading) l/rriQ corrections and three of the four l/mq corrections. 
A similar approach has already been taken for B Dlv 



Our result in the quenched approximation, J^b-^d*{^) = 0.913; 



rO.024+0.017 
-0.017-0.0305 



is consistent 

with results based on other ways of treating non-perturbative QCD. By using the quenched 
approximation we are able to gain control over all other uncertainties. Note, however, 
that the second error bar incorporates (among others) our estimate of the uncertainty from 
quenching. Furthermore, despite the shortcomings of the quenched approximation, it is not 
less rigorous than competing determinations of JFg_^/)*(l), which use either non-relativistic 
quark models or a subjective estimate of the "excitation contribution". With recent mea- 
surements of \Vcb\J^B^D*{X) from CLEO |^ 



the LEP experiments |H9[, and Belle fSO 



our 



result implies 



I • 1 



quark model 



sum rule 



lattice QCD 



0.85 0.90 0.95 



FIG. 6. Comparison of determinations of h^^ (1) = J-b^d* (1 
the non-relativistic quark model a zero-recoil sum rule |4 



with non-perturbative input from 
and quenched lattice QCD. 
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45.9±2.4lli [48] 
lO^\Veb\ = { 38.7 ±1.811:^ PI , (8.1) 



39.3±2.5ii:i g 

where the second, asymmetric error comes from adding all our uncertainties in quadrature. 
Here we have included the QED correction to J-'b^d*{^) of +0.007. 

Since several groups have started partially unquenched lattice calculations of spectrum 
and decay constants, we conclude with some remarks on the prospects for J^b^d*{^)- In 
this context, "partially quenched" means that the valence and sea quarks have different, and 
separately varied, masses. The analysis presented here shows that the double ratios bring 
the statistical precision under control, and that fitting the heavy-quark mass dependence is 
straightforward. Two of our larger systematic uncertainties will improve simply by includ- 
ing dynamical quarks. First, the self-consistent determination of the heavy-quark masses 
and of will improve. At present, we believe the quenching bias in a^, which affects the 
short- distance contribution, to be the largest source of uncertainty from the quenched ap- 
proximation. Second, partially quenched numerical data are enough to extract the physical 
result, because one can use the recently derived result in partially quenched chiral pertur- 



bation theory |51 



The other two main sources of systematic uncertainty are the lattice spacing dependence 
and the matching of lattice gauge theory to HQET and QCD. The former is mostly a 
matter of computing. Indeed, our present estimate may be conservative, as it is driven 
by the coarsest lattice. To decrease the matching uncertainties, one must calculate the 
normalization factor to two loops and calculate the l/mq corrections to one loop. The 
latter is not quite as hard as it might seem. Heavy-quark symmetry protects the needed 
matrix elements, so one only needs the one-loop calculation of the chromomagnetic term 
in the effective Lagrangian (a l/mq term) and the l/mg and mixed \/mjn}j terms in the 
currents. (An alternative to perturbation theory would be to develop a fully non-perturbative 
matching scheme for heavy quarks, including the corrections.) 

With the improvements from unquenched simulations, a more detailed study of lattice 
spacing dependence, and higher order matching calculations, it is conceivable that the error 
on J^B^D'iX) could be brought to or below 1%. At this level, it would become crucial to 
compute, possibly by similar methods, the slope and curvature of Tb^d*{w) near w = 1. 
Then the determination of \Vcb\ would not only become very precise, but also truly model- 
independent. 
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